Interaction of a Nanomagnet with a Weak Superconducting Link 
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We study electromagnetic interaction of a nanomagnet with a weak superconducting link. Equa- 
tions that govern coupled dynamics of the two systems are derived and investigated numerically. 
We show that the presence of a small magnet in the proximity of a weak link may be detected 
through Shapiro-like steps caused by the precession of the magnetic moment. Despite very weak 
magnetic field generated by the weak link, a time-dependent bias voltage applied to the link can 
initiate a non-linear dynamics of the nanomagnet that leads to the reversal of its magnetic moment. 
We also consider quantum problem in which a nanomagnet interacting with a weak link is treated 
as a two-state spin system due to quantum tunneling between spin-up and spin-down states. 

PACS numbers: 75.75.Jn, 74.50.+r, 75.45.+J, 03.67.Lx 



I. INTRODUCTION 

Josephson junctions and magnets in a close prox- 
imity of each other can be coupled through var- 
ious mechanisms. Static properties of supercon- 
ductor/ferromagnet/superconductor (S/F/S) Josephson 
junctions have been intensively studied in the past but 
not as much attention has been paid to the coupled dy- 
namics of the magnetic moment and the tunneling cur- 
rent. The effect of superconductivity on ferromagnetic 
resonance in such junctions has been recently observed by 
Bell et al 1 who attributed their observation to the prox- 
imity effect^. Theory that may be relevant to this exper- 
iment has been worked out by Buzdin who computed the 
phase shift in the Josephson junction arising from the 
Rashba-type spin-orbit coupling 3 - and studied the cou- 
pled dynamics of the magnetization and the Josephson 
current due to this mechanisms. Dynamical proximity ef- 
fect generated by the precession of the magnetization in 
an S/F/S junction has been also investigated by Houzel£ 
Shapiro steps in the FV curve of the S/F/S junction, re- 
lated to the ferromagnetic resonance, have been reported 
by Petkovic et al£ who also provided theoretical argu- 
ments favoring a purely electrodynamic nature of the ef- 
fect in their experiment. 

Coupling of Josephson junctions to individual spins in- 
side the junction have been also intensively studied in the 
past. The theory traces back to the works of Kuliki and 
Bulaevskii et al& who elucidated the effect of spin flips 
on the tunneling current. More recently, Nussinov et 
;iti— . using Keldysh formalism, demonstrated that su- 
perconducting correlations drastically change dynamics 
of a spin inside a Josephson junction. Josephson cur- 
rent through a multilevel quantum dot with spin-orbit 
coupling has been studied by Dell'Anna et alii. Depo- 
sition of a single magnetic molecule in a SQUID loop 
has been attempted^ 2 - and theoretical treatments of the 
Josephson current through such a molecule have been 
propose d 13 ' 14 . Samokhvalov 1 ^ considered formation of 
vortices in a Josephson junction by a magnetic dot. Spin- 
orbit coupling of a single spin to the Josephson junction 



FIG. 1: Color online: Nanomagnet near a weak supercon- 
ducting link. 



has been studied by Padurariu and Nazarov 16 in the con- 
text of superconducting spin qubits. Somewhat related 
to single spins are also studies of two-level systems inside 
Josephson junctions^. 

Early experimental research on coupling of magnetic 
nanoparticles to microSQUIDs has been reviewed by 
Wernsdorfer— who also reviewed recent progress made 
due to the development of nanoSQUIDs 19 . Thirion et 
al2£, demonstrated the possibility of switching of the mag- 
netization of 20nm Co nanoparticles in a dc magnetic 
field by the radio-frequency pulse generated by a mi- 
croSQUID. They measured the angular dependence of 
the switching field and reproduced the Stoner-Wohlfarth 
astroid 21 for a single nanoparticle. Further miniatur- 
ization of such systems has been achieved using carbon 
nanotubes-2 2 - and nanolithography assisted by the atomic 
force microscope 23 . Such systems utilizing single mag- 
netic molecules have been proposed as possible ultimate 
memory units and as elements of quantum computers 24 . 

In this paper we consider a nanomagnet located close 
to a weak link between two superconductors, see Fig. [1] 
The position of the nanomagnet is away from the path of 
the tunneling current, so that the interaction between the 
two systems is considered to be of purely electromagnetic 
origin. The mechanism of the interaction is conceptually 
similar to that argued in the experiment of Rcf. 6. The 
magnetic field of the nanomagnet alters the Josephson 
current flowing through the link, while the magnetic flux 



generated by the Josephson junction acts on the mag- 
netic moment of the nanomagnet. From mathematical 
point of view the dynamics of this problem resembles the 
dynamics studied in Ref. |4|. The differences stem from 
different geometry, different interaction, and finite nor- 
mal resistance of the weak link that we allow within the 
RS J model. The attractiveness of the problem that deals 
with purely electromagnetic interactions is in the absence 
of the unknown parameters. We hope that this will as- 
sist experimentalists in designing Josephson junction - 
nanomagnet systems with desired properties. 

Dynamical equations describing the system depicted 
in Fig. [T] arc derived in the next Section. Small oscilla- 
tions of the magnetic moment of the nanomagnet caused 
by a constant voltage applied to the weak link are stud- 
ied in Section IIIII We show that such oscillations can 
produce Shapiro-like steps in the FV curve without any 
external ac voltage applied to the link. Non-linear dy- 
namics of the nanomagnet due to a voltage pulse applied 
to the link is studied in Section [TV] We show that us- 
ing a specific slow time dependence of the voltage pulse 
one can reverse the magnetic moment of the nanomag- 
net. The remarkable feature of this process is that the 
reversal can be achieved despite the fact that the mag- 
netic field generated by the link is small compared to the 
switching field determined by the magnetic anisotropy. 
This is similar to the effect of the RF field demonstrated 
experimentally in Ref. [2fJ. The reversal occurs due to 
the pumping of spin excitations into the nanomagnet by 
the ac-field of the oscillating tunneling current. The final 
part of the paper studies electromagnetic interaction of 
the weak link with a quantum two-state system formed 
by tunneling of the nanomagnet 's spin between up and 
down orientations. Quantum dynamics of this system 
is derived in Section [V] We show that it provides the 
simplest realization of a Josephson junction - spin qubit 
suggested in Ref. [H (see also Ref. [3 and Ref. |H). Our 
conclusions and suggestions for experiment are summa- 
rized in Section [VTl 



II. THE MODEL 

We consider a system depicted in Fig. [TJ Nanomag- 
net of a fixed-length magnetic moment M is located at 
a distance a from the center of the weak superconduct- 
ing link of length L. The nanomagnet is assumed to be 
rigidly embedded in the solid matrix of the link. In the 
presence of the external magnetic field Bo, the energy of 
the nanomagnet is then given by 



£ M = K(M) - M ■ B : 



(1) 



where K(M.) is the energy of the magnetic anisotropy 
that depends on the orientation of M with respect to the 
body of the magnet. 

Neglecting the capacitance of the weak link, the energy 
of the link can be written as 



where 7 is the gauge invariant phase, Ej = HI c /(2e) is 
the Josephson energy, and I c is the critical current of the 
link. Note that Ej depends on the external field Bo- 
Time derivative of 7, 



dj _ 2eV(t) 

~dl ~ h~' 

is proportional to the total voltage, 



V(t) = J dr-E(r,t), 



(3) 



(4) 



across the link. Here E is the electric field and integration 
goes from one end to the other end of the link. 

For the link biased by the external voltage Vo(t) one 
has 



where 



and 



7 = 7o + 7-4 , 



d 7o = 2eV (t) 
dt h 



7-4 = --— / dr- A(r,t) 



(5) 



(6) 



(7) 



Here $0 = 2irh/(2e) is the flux quantum and A is the 
vector potential. 

In our problem the vector potential A is formed by two 
additive contributions: 



A = A B + A 



Here 



b t -«-M • 



f 



A B = -(B x r) 



(8) 



(9) 



is the vector potential created by the external field and 

^0 M x r 



A M = 



4-7T r 3 



(10) 



is the vector potential created at a point r from the nano- 
magnet assuming that the latter is small compared to all 
other dimensions of the problem. The voltage 



V A = 



h_djA 
2e dt 



(11) 



£,J = -Ej COS 7, 



(2) 



is the electromotive force induced in the link by the time- 
dependent magnetic field generated by the rotating mag- 
netic moment. 

The dynamics of the magnetic moment is given by the 
Landau-Lifshitz equation: 

dM n 

— = 79 MxB e// -^-| 7ff |Mx(MxB e// ) (12) 



where j g is the gyromagnetic ratio for M, r\ is a dimen- 
sionless damping coefficient, and 



Under the above assumptions, substitution of Eq. ([TO 
into Eq. ([7]) gives 



B 



eff 



d£_ 



(13) 



7a = -kM z 



k = 



2-k 



QoaVlX 



(18) 



is the effective field acting on M, with £ being the total 
energy of the system. For £ = £m + £.j one has 



B 



f)T*C F) P^* 

' ff=B °~dM +IcS[ni dMj 1 dr - A M{v,t). (14) 



It is easy to see that the last term in this expression 
equals the magnetic field Bj created by the tunneling 
current I = / c sin7 at the location of the nanomagnet. 
Indeed, substituting into this term Aj\/ of Eq. (TTUJ) and 
rearranging the mixed product of the vectors, one obtains 
for the last term in Eq. (fT4|) 



J c sin7 



Mo d 

4tt<9M 



2 r x dr 

M — 



Mo 

4tt 



2 dl x r' 



= B, 7; 



where r' = — r is the radius- vector pointing from the 
element of the current to the position of the nanomagnet. 
So far we have not considered the normal current 
through the weak link. If the resistance of that link, 
R, is finite, the total current through the link is 



I = I c sin 7 - 



V V 

— = l r sin 7 H 

R c ' R 



ft d-f A 
2eR dt 



(16) 



This expression should replace I c sin 7 in the expression 
for the effective field, so that in the limit of I c — > the 
field given by the last term in Eq. (fT4"|) would be the 
field generated by the normal current In — V/R. Note 
that this field can be formally obtained from Eq. (TlU)) by 
adding the corresponding Zeeman term, 



£z = -In / dr ■ A 



(17) 



to the total energy. 



III. LINEAR APPROXIMATION AND 
SHAPIRO-LIKE STEPS 



Contribution of the weak link to the effective field is 

Bj = kEj r 17 + zr 2^i?^rJ ^ ' (19) 

where e z is the unit vector in the z-direction. Lineariza- 
tion of Eq. (fT^J) then gives the following equation for the 
perturbation of the magnetic moment in the z-direction: 



dm z 



df 2 



■2frr g Bo-^+^B^m 2 



h-f.BvUnE, J I :- ^-i) sin (^L-km z i^n 



73-O0 



Bn dt 



where 



i:\ 2 e lg M 



7] = 71+ [YeJ 2R 



(21) 



is the damping coefficient renormalized by the additional 
channel of dissipation due to normal (eddy) currents gen- 
erated by the rotating magnetic moment. 

According to Eq. (fi"6)) the total current is given by 



I = I,d*{*m-km. 



Vq kh dm z 



R 2eR dt 



(22) 



It has an ac component and the dc component, J, that 
one can obtain by averaging Eq. (|22|) over the oscillations. 
We want to compute the dependence of I on Vb that is 
associated with the I-V curve of the weak link. It is 
convenient to introduce 



2eI c R B Ej 

W ' = V°> ^ = ^T' Eb = T> ^E~ b 



and to switch to dimensionless variables: 

_ M - B _ 

M= , B = — -, t = u a t, 

M S 9 

V = 2 -^, I=L k = kM . 

TlWn I C 



(23) 



(24) 



In this Section we shall assume that deviations of 
the magnetic moment from its equilibrium orientation, 
caused by the interaction with the Josephson junction, 
are small. This will allow us to treat the Landau-Lfshitz 
equation in the linear approximation. For certainty, we 
choose the external magnetic field Bo and the equilib- 
rium magnetic moment Mo in the direction parallel to 
the line connecting the leads 1 and 2, which is the y- 
direction in Fig. [TJ To make the problem more tractable 
we shall also assume in this Section that the applied field 
is large compared to the effective field due to magnetic 
anisotropy, so that the latter can be neglected. 



Note that uj g is the precession frequency for the magnetic 
moment in the absence of interaction with the supercon- 
ducting link. 

In terms of the above variables the equations for I and 
m z become 



I = sin(y t - krh z ) + ^-(v - k^ 
uj r \ dt 



(25) 



i 2 m z dm z ( d , 



(26) 
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FIG. 2: Color online: Dependence of I on Vq in the presence 
of the nanomagnet (purple) and without the magnet (blue) for 
e = 0.1, k — 0.1, fj fa r\ — 10 -4 , lu 3 /lor <C 1. The inset shows 
the difference between the two curves. 



In these equations the dimensionless parameter k can be 
small or large, depending on the size and the location of 
the magnet. The ratio uj g /u)R can also be small or large 
depending on the resistance R. The parameter e roughly 
equals the ratio of the field created by the critical current 
at the location of the nanomagnet and the external field. 
In practical situations this ratio will always be small, 
thus, justifying the linear approximation for rh z away 
from resonance, Vq — 1, and at the resonance for not 
very small rj. In the case of a very narrow resonance (very 
small fj) one should employ the non-linear approximation 
based upon the full Landau-Lifshitz equation. 

The dependence of (/) on Vq, computed numerically, 
is shown in Fig. [2] Shapiro-like steps at Vq = 1 and 
Vq = 2 are apparent. They appear due to same physics 
as the conventional Shapiro steps, with the field of the 
precessing magnet playing the role of the rf field. The 
half-Shapiro step that can be seen at Vq = 0.5 appears 
when one solves the full Landau-Lifshitz equation 
instead of the linearized equation. Fig. [2] illustrates the 
principal possibility to detect the presence of a small 
magnet in the vicinity of the weak link by measuring its 
I-V curve. 



(with K being a constant and V being the volume of the 
magnet) in a zero external field. The effective field from 
this term in the energy is 



B 



(A) 
eff 



KMy 

V 



(28) 



The problem of coupling of the weak link to small oscil- 
lations of M around M directed along the anisotropy 
axis becomes identical to the problem studied in the pre- 
vious section if one replaces Bq with KMq/V. If both 
are present Bq in the above formulas should be replaced 
with Bq + KMq/V. 

The dc magnetic field that would be required to switch 
the magnetic moment to the opposite orientation along 
the anisotropy axis is Bq = KMq/V. In all practical 
situations the magnitude of the oscillating magnetic field 
produced by the weak link at the location of the mag- 
net will be hopelessly small compared to that field. The 
question, however, arises whether the ac field produced 
by the oscillations of the Josephson current can pump 
spin excitations into the magnet at a rate sufficient to 
reverse its magnetization. We shall see that this may, 
indeed, be practicable, thus invoking the possibility of a 
magnetic memory unit operated by voltage pulses. 

When the effective field is dominated by the magnetic 
anisotropy the roles of the parameters to g and Eb are 
played by 



JgKMp 

V 



Er 



KMq 
kV 



(29) 



To simplify our formulas we consider in this Section the 
limit of a very large normal resistance, so that we can 
neglect the normal current through the link. (This as- 
sumption is unessential for our conclusions, though, and 
the calculation can easily be generalized to the case when 
the normal current is present.) Under this assumption 
the non-linear dynamics of the magnetic moment is de- 
scribed by the dimensionless Landau-Lifshitz equation 



dM 



= M X B e // - 7?M X (M X B e //) (30) 



IV. NON-LINEAR DYNAMICS AND 
MAGNETIZATION REVERSAL 



with dimensionless 



B e // = M y e y +esin(Vo* - kM z )e z 



(31) 



In this Section we demonstrate the possibility of a re- 
versal of the magnetic moment of the nanomagnet by 
using a specific time dependence of the bias voltage ap- 
plied to the weak link. This problem involves a non- 
linear dynamics described by the full Landau-Lifshitz 
equation. Consider a nanomagnet with uniaxial mag- 
netic anisotropy 



K(M) 



2V v 



(27) 



Numerical solution of Eq. (|30|) for a time-linear voltage 
pulse is shown in Fig. [3] As the effective field decreases 
in the course of the reversal, the ac field generated by 
the oscillating tunneling current continuously pumps spin 
excitations into the nanomagnet. This leads to the full 
reversal of the magnetic moment. We find numerically 
that the reversal only occurs at e > r\. Another observa- 
tion is that the time needed for the reversal is inversely 
proportional to e. Practical implications of these findings 
are discussed in Sec. I VII 




FIG. 3: Color online: Reversal of the magnetic moment of the 
nanomagnet by linearly decreasing bias voltage V = 1.5 — 
0.00075£. The values of the parameters are e = k — 0.05, 
77 = 0.01. 



V. NANOMAGNET AS A TWO-STATE 
QUANTUM SYSTEM 

In this Section we will treat nanomagnet as a fixed- 
length quantum spin S, rigidly embedded in a solid ma- 
trix. Magnetic anisotropy energy K(M) should now be 
replaced by a crystal-field Hamiltonian. The general form 
of such a Hamiltonian that corresponds to a strong easy- 
axis magnetic anisotropy is 



He 



H l{ +H A 



(32) 



where Hu commutes with S z and H± is a perturbation 
that does not commute with S z . Presence of the mag- 
netic anisotropy axis means that the | ± S) eigenstates 
of S z are degenerate ground states of Ht\ . Operator H± 
slightly perturbs the | ± S) states, adding to them small 
contributions of other \ms) states. We shall call these de- 
generate normalized perturbed states iV'is)- Physically 
they describe the magnetic moment of the nanomagnet 
looking in one of the two directions along the anisotropy 
axis. Full perturbation theory with account of the de- 
generacy of Hs provides quantum tunneling between the 
\t/j±s) states^. The ground state and the first excited 
state are even and odd combinations of \ip±s) respec- 
tively, 



They satisfy 



with 



*T = ^= (Ilk?) ± |tf-S» 



H S ^± = E±if>± 



E-l- — E- 



A 



(33) 



(34) 



(35) 



being the tunnel splitting. The latter is typically small 
compared to the distance to other spin energy levels, 
making the two-state approximation rather accurate at 
low energies. For, e.g., biaxial magnetic anisotropy, 



H s = -DS 2 Z + dSl with d < D, the splitting of the 
lowest energy level appears in the iS-order on d/ D, while 
the distance to the next level equals (25 — 1)D. 

Since the two low-energy spin states of quantum nano- 
magnet are superpositions of \ip±s), it is convenient to 
describe such a two-state system by a pseudospin 1/2. 
Components of the corresponding Pauli operator <x are 

&x = \lp-s){lps\ + \lps){i>-s\ 

(T y = i\ip-8){ips\-i\'4>s){'<P-s\ 

o z = \$s)l$s\-\1>-s)l$-s\. (36) 

The projection of any operator As — A(S) onto \4>±s) 
states is 

A a = J2 HA s \n)\m)(n\. (37) 

Expressing |i/>±s) via *± according to Eq. (|3^)l . it is easy 
to see from Eq. ([M]) that 

(il>±s\6 S \il>±s) = 0, (il>-s\6 S \1>s) = -A/2 . (38) 

With the help of these relations one obtains from Eq. 
@7J) 



H^=-(A/2)* X 



(39) 



Quantum generalization of Eq. ([2]) with account of Eqs. 
© and (HU) is 

Hj = -Ej cos(ujjt - k/J, B S z ) ■ (40) 

Equations (|3"rj)) and (13"T|) then give 

H^ l) = -Ej[cos{ujjt) + kfi B Sa z sin(wjt)] , (41) 

where 



wj(t) = yVo(t). 



(42) 



The total Hamiltonian of our two-level system is 

H a = -b e// ■ a , (43) 

where 

b e// = (A/2)e x + kEj sin(LOjt)e z . (44) 

Here, as before, k = kM = k^sS, so that the energy 
kEj roughly represents the strength of the interaction 
of the magnetic flux of the junction with the spin S. In 
practice this interaction can be greater or smaller than 
A. 

The above system represents a simple realization of 
the spin qubit proposed by Nazarov^. Quantum states 
of such a qubit are described by the wave function 



* = C s \ips) + C_s|^_ s ) 



(45) 
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FIG. 4: Color 
C S (0) = 1. 



online: Rabi oscillations of \Cs(t)\ 2 at 2kEj/A — 0.1 for two constant voltages with the initial condition 



The Schrodinger equation for |\&) is 



d* 

in—— 
dt 



(46) 



Here |C±s| 2 is the probability for the spin to look up or 
down along the z-direction, with |C + s| 2 
Introducing 



|C_ 



= 1. 



C±s(t) = C ±s (t)exp 



ikEj 
2ft 



dt'sin[uj(t')t' 



we obtain from Eq. (|4l)|) 

lh dt°- S = ^° S 
ih—Cs = — kEjCs sm(ujt) 



(47) 






(48) 



In terms of dimensionless time, t — At/ (2ft), the resulting 
equation for C# is 



dP 



Cs 



.2kE j d 

A ~dl 



sin(^!r, 



+ C s = (49) 



The effect of the bias voltage becomes especially pro- 
nounced at Vo satisfying ujj = (m/n)A, where m and 
n are integers. Fig. [4] shows Rabi oscillations of the 
probability to remain in the initial spin-up state for two 
different bias voltages each satisfying one of the above 
conditions. 



VI. DISCUSSION 

We have studied electromagnetic interaction between 
a weak superconducting link and a small magnet placed 
in the vicinity of the link. Three problems have been 
considered: Shapiro-like steps in the I-V curve generated 
by the magnet, the reversal of the magnetic moment by 



a time-dependent bias voltage, and Rabi oscillations of 
the quantum spin induced by a constant voltage. 

For the first two problems the strength of the inter- 
action is determined by the parameter e = Ej/Eb- By 
order of magnitude e represents the ratio of the magnetic 
field generated by the tunneling current and the effec- 
tive field, B e ff, acting on the magnetic moment due to 
magnetic anisotropy and the applied external field. We 
show that precession of the magnetic moment generates 
Shapiro-like steps in the I-V curve of the superconduct- 
ing weak link. The possibility to observe the first Shapiro 
step at Vq = 1 (and also the peak at Vq = 0.5 due to non- 
linearity) appears quite realistic. Note that the first step 
scales down linearly with e when decreasing e, the second 
step at Vq — 2 scales as e 2 , and so on. Thus, for small e, 
higher steps may be more difficult to see in experiment. 

A remarkable observation is that despite the weakness 
of the field generated by the tunneling current of the link, 
for a certain time dependence of the bias voltage it can 
effectively pump spin excitations into the magnet, lead- 
ing to the reversal of its magnetic moment. The damping 
constant r) — 0.01 was chosen for simulations of the rever- 
sal. This value is realistic for magnetic nanoparticles 27 . 
We find that condition e > r\ is required for the reversal, 
which must be important for experiment. The parame- 
ter e determines the number of cycles in the precession 
of the magnetic moment that leads to the reversal of 
the moment. In our numerical simulations that number 
roughly scaled as 1/e. For e = 0.05 used to obtain the 
plot shown in Fig. [3] the time required to reverse the mo- 
ment was close to lO 3 ^!" 1 . For, e.g., uj g ~ 10 11 s — x this 
would provide a reversal in ten nanoseconds. The linear 
time dependence of the bias voltage in Fig. [3] was cho- 
sen to maintain the condition of continuous pumping of 
spin excitations into the magnet. Smaller e would require 
slower time dependence of Vq , which should not be diffi- 
cult to satisfy in experiment. However, smaller e would 
require smaller r\ due to the condition e > r\. Also, the 
smaller is e the more sensitive the time evolution of the 



magnetic moment becomes to the time dependence of the 
voltage. A slight change in that dependence is sufficient 
for the moment to bounce back to the initial direction 
after reaching the top of the anisotropy barrier. 

In the quantum problem, the parameter e is no longer 
relevant. The relevant parameter becomes the ratio of the 
Zeeman interaction of the spin with the field of the tun- 
neling current and the tunnel splitting A. This parame- 
ter can be small or large depending on the splitting. Rabi 
oscillations of the spin are strongly affected by the bias 
voltage. The most noticeable effect appears at V satisfy- 
ing one of the resonant conditions eVo = (m/n)A, where 
m and n are integers. At such resonances the behavior 
of the probability to find the spin in up or down config- 



urations is very different from the off-resonance behav- 
ior. This demonstrates the principal possibility to elec- 
tromagnetically manipulate a nanomagnet - weak link 
qubit by the voltage applied to the link. 
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